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We propose a scheme of multilayer thermoelectric engine where one electric current is coupled 
to two temperature gradients in three-terminal geometry. This is realized by resonant tunneling 
through quantum dots embedded in two thermal and electrical resisting polymer matrix layers 
between highly conducting semiconductor layers. There are two thermoelectric effects, one of which 
is pertaining to inelastic transport processes (if energies of quantum dots in the two layers are 
different) while the other exists also for elastic transport processes. These two correspond to the 
transverse and longitudinal thermoelectric effects respectively and are associated with different 
temperature gradients. We show that cooperation between the two thermoelectric effects leads to 
markedly improved figure of merit and power factor which is confirmed by numerical calculation 
using material parameters. Such enhancement is robust against phonon heat conduction and energy 
level broadening. Therefore we demonstrated cooperative effect as an additional way to effectively 
improve performance of thermoelectrics in three-terminal geometry. 

PACS numbers: 73.63.-b,85.80.Fi,85.35.-p,84.60.Rb 


I. INTRODUCTION 

Harvesting usable energy from wasted heat using ther¬ 
moelectrics has been attracting a lot of research interest 
Much efforts have been devoted to bulk materials, mak¬ 
ing them mature thermoelectric systems that are al¬ 
ready useful in industrial technologiesRecently there 
is a trend of incorporating nanostructures to further im¬ 
prove the performance of thermoelectric materialsi^^— 
Successful examples have been achieved in many mate¬ 
rials/structures along this directional^— Theoretical stud¬ 
ies have demonstrated that the thermoelectric properties 
of individual nanostructures can be much better than 
the bulk— Experimental efforts have pushed forward 
the measurements of thermoelectric properties of indi¬ 
vidual nanostructures There are also studies try¬ 
ing to fill the gaps between the thermoelectric proper¬ 
ties of individual nanostructures and their assemblies'^— 
as well as attempts to improve thermoelectric perfor¬ 
mance by tuning the shape and organization patterns 
of nanostructures—^ 

In addition to material and structural aspects, geome¬ 
try also plays an important role in thermoelectric appli¬ 
cations. For example, transverse thermoelectricsi^ take 
advantages of accumulating temperature difference in one 
direction while generating electric current in the perpen¬ 
dicular direction. Geometric separation of the electric 
and heat flows facilitates special functions. For example, 
thermoelectric cooling and engine can be realized using a 
single type of carrier doping (i.e., without serial connec¬ 
tion between n- and p-type thermoelectric components) 
via transverse thermoelectric effect—^ Recently a related, 
but different, thermoelectric effect is found in meso¬ 
scopic thermoelectrics in three-terminal geometry— 
Researches in this direction is pioneered by the theory 


of Edwards et al^ and the later experiments— ^ The un¬ 
derlying physics is illustrated in Fig. [T] (see also Ref. [l^: 
excess population of phonons can induce an electric cur¬ 
rent during inelastic transport processes. Heat and elec¬ 
tric flows are geometrically separated since heat is car¬ 
ried by the phonons flowing from/into the phonon bath. 
This picture can be generalized to inelastic transport pro¬ 
cesses assisted by other elementary excitations, such as 
photons^S, electron-hole excitation s^^d^ and magnonsi^ 
Besides the quantum dots (QDs) can be replaced with 
any conductors given that the carrier energies at the left 
and right conductors are considerably different, which 
can be realized by two low-dimensional structures (e.g., 
quantum wells^^ or wires), or a barrier^ or a band gapi^ 

Microscopic analysi a^^d'^~^^d^~ iS^ indicates that the 
performance of each individual nano-scale three-terminal 
thermoelectric (3T-TE) device is promising. Experi¬ 
ments have demonstrated the effectiveness of 3T-TE 
cooling at submicron scalei^ In this work we focus on 3T- 
TE systems based on the structure illustrated in Fig. 2. 
This structure was initially proposed by Edwards et al^ 
and later explored experimentally in Ref. HO for cool¬ 
ing of electrons at cryogenic temperature. Recently Jor¬ 
dan et al. extend the idea to thermoelectric engine with 
layered self-assembled QDs where considerable electrical 
current density could be obtained due to contributions 
from many parallel quantum tunneling channels^. This 
proposal significantly improve the potential for thermo¬ 
electric energy harvesting of the original idea (An ex¬ 
tended idea of replacing the QD layers by quantum wells 
is presented in Ref. [28l) . Here we exploit the same struc¬ 
ture of Jordan et al. to study cooperative effects between 
the longitudinal and transverse thermoelectric powers. 

In Fig. Ha) the electronic cavity, as well as the source 
and the drain, are highly conducting layers made of 
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grown on the surface of semiconductors;^ or as we pro¬ 
posed here, core-shell QDs embedded in (undoped) poly¬ 
mers with low thermal and electrical conductance42r— 
In such an assembly scheme the total electric current is 
the sum of the electric currents in each nano-scale 3T- 
TE device (i.e., a pair of QDs). High power density can 
be prompted by high density of QDs which can reach 
^ 10^^ cm“^ for a single layer (about 1 nm thick)i^ 
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FIG. 1. (Color online) Illustration of a three-terminal ther¬ 
moelectric system (stndied in Ref. [l^. (a) An electron first 
tunnels from the source to the quantum dot (represented by a 
red dot) with energy Ei and then hops to the other quantum 
dot with energy E 2 by absorbing a phonon from the phonon 
bath. After that the electron tunnels into the drain elec¬ 
trode. The tunneling processes are elastic, while the hoping 
between quantum dots must be assisted by a phonon with en¬ 
ergy E 2 — El > 0. The temperatures of the source, drain and 
phonon bath are Ts, Td, and Tp, while the electrochemical 
potentials of the source and drain are fis and fj,d, respectively. 
When the phonon bath has the highest temperature among 
the three reservoirs, excess number of phonons prefer phonon 
absorption processes and hence induces an electric current 
from the source to the drain, (b) In realistic situations the 
phonon bath is connected with the (insulating) substrate sup¬ 
porting the quantum dots while the two electrodes are sus¬ 
pended to be isolated from the phonon bath. 



heavily-doped semiconductors (e.g., heavily-doped sili¬ 
con or GaAs). In between those layers, there are two 
highly resisting layers with high thermal and electrical 
resistance. Each of the resisting layer is embedded with a 
QD through which electrons can tunnel between the cav¬ 
ity and the electrodes. The resonant tunneling through 
QDs are responsible for the transport. The energy lev¬ 
els in the left and right QDs are Ei and E 2 , respectively. 
When El ^ E 2 , an electron transmitting from the source 
to the drain takes a finite amount of energy, E 2 — Ei, 
from the cavity. To reach steady states the cavity must 
exchange energy with the phonon bath [see Fig. [ll[b)]. 
The transverse thermoelectric effect is manifested as the 
fact that an electric current drives a heat current from 
the phonon bath [see Fig. [21[b)], and vice versa. 

Following Ref. [13, the scheme to assemble the nano¬ 
scale devices into a macroscopic device is straightforward: 
2D arrays of QDs can be placed in the resisting layers [see 
Fig. me)]. This can be realized by self-assembled QDs 


FIG. 2. (Golor online) Illustration of the 3T-TE device stud¬ 
ied in this work, (a) Schematic of a possible realistic set-up. 
The layered structure consists of three highly conducting re¬ 
gions: the source, cavity, and drain. In between them there 
are two highly resisting layers (yellow regions) with very low 
thermal and electrical conductance. A QD is embedded into 
each resisting layer to allow electrical and thermal conduc¬ 
tion through it. The cavity is connected with a (insulating) 
phonon bath. The connection allows efficient energy exchange 
between them, (b) Illustration of the microscopic processes. 
Garrier distributions in the source and the drain are deter¬ 
mined by the chemical potentials and temperatures. Elec¬ 
trons transmit through QDs via resonant tunneling. When 
energies of the two QDs, Ei and E 2 , are different, energy 
exchange between the cavity and the phonon bath is neces¬ 
sary for the establishment of the steady state transport, (c) 
A scheme for building macroscopic devices where many QDs, 
forming 2D arrays, are embedded in the resisting layers. 

To date all works on 3T-TE systems focus on exploiting 
the transverse thermoelectric effect. However, there is 
also a longitudinal thermoelectric effect in the system: 
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the temperature difference between the two electrodes 
can also induce an electric current. A full description of 
thermoelectric transport in 3T-TE systems is given by 
the phenomenolo gica l equation (similar equations were 
found in Refs, fm Il9l. 1^. and [2^ 

f I \ / G Li L2 \ / V \ 

/qi = Li K, Ku ATi/T , (1) 

\IQ 2 J V L 2 Ku K 2 J \ AT2/T j 

where Li and L 2 represent the longitudinal and trans¬ 
verse thermoelectric effects, respectively. Iqi and Iq 2 
stand for the heat currents leaving the source and the 
phonon bath, respectively. The two temperature differ¬ 
ences are ATi = Tg — Tj^ and AT 2 =Tp — T^ with T^, 
and Tp being the temperatures of the source, the drain, 
and the phonon bath, respectively. 

In this work we show that, if both the longitudinal 
and transverse thermoelectric effects are exploited simul¬ 
taneously, due to cooperation between the two, the ef¬ 
ficiency and power can be considerably improved. The 
cooperative effect originates deeply from the nature of 
three-terminal thermoelectric systems: the two thermo¬ 
electric effects are correlated with each other, or in other 
words, the electrical current are simultaneously induced 
by the two different temperature gradients. Simplified 
geometric interpretation is that the electric currents in¬ 
duced by the two thermoelectric effects can be parallel or 
anti-parallel. In the former case the two effects add up 
constructively, leading to enhanced thermopower and ef¬ 
ficiency. The phenomenon reflects the cooperative effect 
of two (or more) correlated thermoelectric effects which 
is referred to as the “cooperative thermoelectric effect”. 

The cooperative thermoelectric effect is manifested in 
the proposed thermoelectric device. Using material pa¬ 
rameters we calculate the thermoelectric transport coef¬ 
ficients as well as the figure of merit and power factor. It 
is found that both the figure of merit and the power fac¬ 
tor are considerably improved by the cooperative ther¬ 
moelectric effect. This enhancement is as effective for 
good thermoelectrics as that for bad thermoelectrics. In 
calculation we show that the enhancement induced co¬ 
operative effect changes only slightly when phonon heat 
conductivity or QD energy broadening is increased sig¬ 
nificantly. These results demonstrate that cooperative 
effect is an alternative way to improve the performance 
of thermoelectrics in three-terminal geometry effectively. 

This paper is organized as follows: In Sec. II we es¬ 
tablish thermoelectric transport of the 3T-TE systems 
from microscopic theory. In Sec. Ill we demonstrate the 
cooperative thermoelectric effect in a geometric way. In 
Sec. IV we calculate the thermoelectric transport coeffi¬ 
cients as well as the figures of merit and power factors for 
the longitudinal, transverse, and cooperative thermoelec¬ 
tric effects for the 3T-TE systems using material param¬ 
eters. We conclude in Sec. V. Studies in this work are 
focused on linear-response steady state transport. In¬ 
teresting nonlinear effects^^ could be discussed in future 
works. 


II. MICROSCOPIC THEORY OF 
THERMOELECTRIC TRANSPORT 

In this section we develop a microscopic theory of 
3T-TE transport, following the formalism of Jordan et 
The system is a layered structure with thickness 
^tot — Tg-f Lg[-\-L q-\- ‘ZLqfi where Lsi and Lqd are 

the thickness of the source, the drain, the cavity, and the 
resisting layers, respectively. For realistic design, Lg, L^, 
and Lc is about one hundred nanometers (nm), while Lqd 
is on the order of ten nm. The size of QDs is a few nm. 
The level spacing of the QDs, typically on the order of 
100 meV for core-shell QDs (see Ref. [33), is much larger 
than ksT. We suggest to fabricate serially connected 
unit structures [the structure in Fig. UKa)] to scale the 
device up to a fully three-dimensional macroscopic de¬ 
vice which can be implemented via layer-by-layer growth 
methods. 

The Hamiltonian of the system is written as 


H = Hg + Hd + H, + Hqd + mr,u (2) 

where i/g, Hd, He, and Hqd are the Hamiltonian of 
the source, the drain, the cavity, and the QD, respec¬ 
tively. Ha = o’ where a = s,d,c denotes the 

k2i2 

source, the drain, and the cavity, respectively, eg = 
with m* being the effective mass of the charge carrier. 

Hqd = 4 + 4 ’ ( 3 ) 

* j 


where the index i and j numerate the QDs in the left 
(£) and right (r) resisting layers, respectively. Hint de¬ 
scribes hybridization of the QD states with the states in 
the source, drain, and cavity. 


Hjn 


= EE^. 

OL — S,C 


d, 

k k.oc 


y yy^ct dj+R.c. 

< ^ ^ j,Oi,k k,ot 


y.—d,( 


k,j 


( 4 ) 

The electric and thermal currents through the left resist¬ 
ing layer from the source to the cavity are given by^ 


Ie,i =‘^JdE ti{E)[fg{E) - fe{E)], (5a) 

/q,, = ^J dE(E - ME)[fs{E) - fe{E)], (5b) 

respectively. The factor of two comes from the spin de¬ 
generacy of the carriers, e is the electronic charge, /g 
and fc are the carrier distribution functions of the source 
and the cavity, respectively. They are determined by the 
temperatures of the source Tg and the cavity Tc as well 
as by their electrochemical potentials g(g and /ic- Note 
that because the voltage and temperature gradients are 
mainly distributed at the two resisting layers, as a good 
approximation, one can assign uniform chemical poten¬ 
tials and temperatures to the source, drain, and cavity re¬ 
gions. The energy dependent transmission through QDs 
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is given tt{E) = ti{E) with 

, , fi^^s4E)ri^l{E) 

U{E) = - -2 > (6) 

{E-Ef'>Y + !^{Ts,^[E)+T^l{E)) 

where 

k 

= (7b) 

k 

are the energy-dependent tunneling rates from the QD i 
to the source and the cavity, respectively. The electric 
and thermal currents from the drain to the cavity can be 
obtained by the replacements: z —>■ j, £ —r, and s ^ d. 
Introducing 

Gi = ^ J dE U{E)f{E)[l - f{E)\ (8a) 

L, = J dE{E - ^^)U{E)f{E)[l - f{E)] (8b) 

K, = J dE{E - ^ifu{E)f{E)[l - f{E)\ (8c) 

with f{E) being the equilibrium carrier distribution, in 
the linear-response regime one can rewrite Eq. as 

le^t = Geifis - Mc)/e + Le{Ts - Tc)/T, (9a) 

lQ,i = - fj.c)le + Ki{Ts - Tc)/T. (9b) 

Expressions for the currents from the drain to the cav¬ 
ity Ie,r and Ig^r can be obtained from the above by the 
replacements i ^ r and s —>■ d. 

Inelastic scatterings, such as the electron-phonon and 
electron-electron scatterings, are crucial for the estab¬ 
lishment of steady states in the cavity. In the concerned 
temperature range, 300 ^ 500 K, those scatterings are 
quite efficient. The heat transfer between the phonon 
bath and the cavity can be made efficient by using ma¬ 
terials with high thermal conductivity to connect them. 
Interface thermal resistance can be reduced if the cavity 
and the phonon bath are made of the same material. We 
assume the thermal conduction between the phonon bath 
and the cavity is efficient and the temperature gradient 
across them and within the cavity is considerably smaller 
than that across the two polymer layers. In this way the 
temperature of the cavity is very close to that of the 
phonon bath and one can approximate that Tc = 

Energy conservation gives -I- Iq^r + lQ,t + IV = 0, 
where is the heat current from the phonon bath to 
the cavity. Therefore there are only two independent 
heat currents In Eq. o the two independent heat 
currents are chosen as Iqi = and Jq 2 = Iq^t- Charge 
conservation, + Ie,r = 0, determines the electrochem¬ 
ical potential of the cavity 

He = {Gi + Gr) ^^GiHe + GrHd + sT ^ {Tg — Tp) 
+L,(Td-Tp)]}. (10) 


Inserting the above into Eqs. m and (jH) we obtain 

,Li 


G = {GJ^ + G-^) 


iN-l 


Li = G- 


L2 = G{^ - ^), 

Ll{Li + Lr) 

Ei2 = — t ;— 7 “^;- Ei, 


Gi 

Ki = Ki- 


Li 


Gi + Gr 


Gi + Gr 
K2 = Ki + Kr- 


{Li + 


(lla) 

(llb) 

(llc) 
(lid) 


Gi + Gr 

To understand these results, we rewrite the transport co¬ 
efficients in terms of average electronic energies, following 
Mahan and Sofo,— 

Li = e-^G{E- H)i, L 2 = e-^G{{E)^-{E)^), (12a) 


7^2 


G 


Ko = 


G 


Kti, 

■ Ktl + Ktr 


Ij'] JjO 

Ki2 = —- Ktl, 


G 


where 


Ktl = e ^Gi 


Ktr — ' 


^Gr 


(12b) 

(12c) 

(13a) 

(13b) 

(14) 


{{E-hY)^-{E-h)] 
{{E-hY)^-{E-h) 

The average in the above is defined as 
,jpn. _ jdEE-Gp{E) 

^ JdEGpiE) 

with G^iE) = ^tp{E)f{E)[l - f{E)], for ^ = i,r 
and n = 0,1, 2. Gi = f dEGi(E) and Gr = f dEGr{E). 
One readily notices from Eq. (USD that Ktl and Ktr must 
be non-negative. 

Eor a macroscopic system with area A the electrical 
conductivity is cr = Glu/A with = 2Zgd + 4 + (G + 0)/2 
being the thickness of an unit structure [the structure in 
Eig.j^Ka)]. Similarly the thermal conductivities are zti = 
KiL/{AT), K 2 = KiluliAT), and «:i 2 = Ki 2 lut{AT). 
The longitudinal and transverse thermopowers are 

(15) 


Si s 


El {E — H)i 


S2 = 


TG eT ’ TG eT 

S 2 is proportional to the energy difference, reflecting that 
it is associated with the inelastic processes. In contrary 
Si remains finite when inelastic processes vanish. 

The total entropy production of the system in the lin¬ 
ear response regime is written as 


dt T ' 


Iqi 


ATi 


Iq2 


An 


T 


(16) 


The second law of thermodynamics, ^ > 0, requires 
that^ 


GKi > Ll, GK 2 >Ll, KiK 2> Kf 


(17) 


as well as that the determinant of the 3x3 transport ma¬ 
trix in Eq. o to be non-negative. Those requirements 
are satisfied for the transport coefficients in Eq. (fl^ be¬ 
cause Ktl, Ktr > 0. 
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III. COOPERATIVE EFFECT: A GEOMETRIC 
INTERPRETATION 


We parametrize the two temperature differences as 


ATi = ATcos 6 l, AT 2 = ATsme. (18) 

The exergy efficiency (or the “second-law efficiency”, see 
Refs. 113 and 113) of the thermoelectric engine is^“— 


_-(w_ yzT + 1-1 

/qi ATi/T + IQ 2 AT 2 IT - V^T+1 + 1 ■ 


The exergy efficiency (“second-law efficiency”) is defined 
by the output free energy divided by the input free 
energy^^— It has been widely used in the studies of en¬ 
ergy conversion in chemical and biological systems since 
its invention about 60 years ago4i Recently it was ap¬ 
plied to thermoelectric systems^. According to Ref. 
the rate of variation of free energy associated with a cur¬ 
rent is given by the product of the current and its conju¬ 
gated thermodynamic force. Hence the denominator of 
the above equation consists of heat currents multiplied 
by temperature differences. It has been shown in Ref. |13 
that the relation between the efficiency of 77 = W/Q 
for heat engine (or 77 = QjW for refrigerator) and the 
second-law efficiency (f is that (f = rj/rjc where rjc is the 
Carnot efficiency. Thus Ioffe’s figure of merit is also ob¬ 
tained starting from the second-law efficiency. At given 
0 the figure of merit is 


ZT = 




Keff - O'S'gfjT 


(19) 


is the figure of merit. Here S'es = Si cos 6 + S 2 sin0 and 
Keff = (A'l cos^ 9 + 2Ki2 sin 9 cos 9 -I- K 2 sin^ 9)1^/{AT). 
Upon optimizing the output power of the thermoelectric 
engine, one obtains^ 

Wmax = ^PiAT)^, ( 20 ) 


and TT < 9 < 37 r/ 2 , ZT is greater than both ZiT and 
ZtT. To understand the underlying physics, we decom¬ 
pose the electric current into three parts I = Iq + Ii +12 
with lo = GV, h = LiATi/T, and h = L 2 AT 2 /T. The 
two thermoelectric effects add up constructively when 
Ii and I 2 have the same sign which takes place when 
0 < 0 < 7 r /2 and tt < 9 < 37 r/ 2 . Fig. El^b) shows the 
power factor versus the angle 9. The power factor is also 
larger when the two currents Ii and I 2 are in the same 
direction. The cooperation of the two thermoelectric ef¬ 
fects thus leads to enhanced figure of merit and output 
power. 



with the power factor 


P = aS^ 


eff- 


( 21 ) 


When 0 = 0 or TT, Eqs. (US and m give the well-known 
figure of merit and power factor for the longitudinal ther¬ 
moelectric effectii^ 


ZiT = 


aSfT 

Ki-aSfT’ 


Pi = a Si 


( 22 ) 


The transverse thermoelectric figure of merit and power 
factor, i.e., 9 = tx 12 or 37 r/ 2 , are given 


ZtT = 


<ySlT 

K2-aSjT’ 


Pt = crSl 


(23) 


Fig.EKa) shows ZT versus the angle 0 in a polar plot for 
a specific set of transport coefficients satisfying the ther¬ 
modynamic bounds in (flTl) . Remarkably for 0 < 0 < 7 r /2 


FIG. 3. (Color online) Polar plot of figure of merit ZT (a) and 
power factor P (b) [in arbitrary unit (a.u.)] versus angle 9. 
At 0 = 0° or 180° ZT and P recover the values for the longi¬ 
tudinal thermoelectric effect (red dots), while at 0 = 90° and 
270° they go back to those of the transverse thermoelectric ef¬ 
fect (green squares). The arrows in the I, II, III, IV quadrants 
label the direction of the currents 7i = LiATi/T (red arrows) 
and I 2 = L 2 AT 2 /T (green arrows). The transport coefficients 
are: L\ — L 2 = 2GkBT/e, Ki = K 2 = 12G(fesr/e)^, and 
7^12 =0. 7o = GV (not shown) is reversed in III and IV 
quadrants for the operation of the thermoelectric engine. 


One can maximize the figure of merit by tuning the 
angle 9. This is achieved at 


de{ZT) = 0. (24) 


We find that the figure of merit is maximized at 0 = 9m 
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with 


tan(0M) = 


L2K1 — L1K12 
L1K2 — L2K12 


S2K1 — S1K12 
S'lAt 2 ~ S2K12 


(25) 


After some algebraic calculation the maximum figure of 
merit is found to be 


ZmT 


G{KiK 2 - Kl^) 
Dm 


(26) 


where Dai = GKiK2-GK^2-K2Ll + 2LiKi2L2-KiLl 

denotes the determinant of the 3x3 transport matrix in 
Eq. ([T]). ZmT is greater than both ZiT and ZtT, un¬ 
less the denominator or the numerator in Eq. (|25l) van¬ 
ishes. Nevertheless it is guaranteed by Eqs. (m and m 
that both the numerator and denominator in Eq. (1251) is 
nonzero when the broadening of the quantum dot energy 
is finite. 

One can also tune 0 to maximize the power factor P 
which is achieved at 


deP = 0 . 

(27) 

The power factor is maximized at 9 = 
9m ^ 9 m) with 

9m (in general 

, tn \ _ T 2 S 2 

tan(G'j72) — — . 

Li bi 

(28) 

The maximum power factor 


Pm = a{Sj + Si) 

(29) 


is greater than both Pi and Pj unless or S 2 is zero. 
If 9m is close to 6m, both the figure of merit and the 
power factor can be improved by the cooperative effect 
simultaneously in certain range of 9. 

The cooperative thermoelectric effect becomes partic¬ 
ularly simple and vivid when K 12 = 0 , Ki = K 2 = AT, 
Li = L 2 = L, and S = L/{TG). In this special case the 
transport equation becomes 


= ( 30 ) 

\Iq2 J \l 0 K J \ AT 2 /T J 


We shall use the following combinations of temperature 
differences 


AT„. = 


ATi AT 2 

7i 


An = 


ATi - AT 2 

7 i 


( 31 ) 


The heat currents conjugate to the above forces are 


f Iqi + ^Q 2 T 

“ V2 • ^ 

The transport equation then becomes 

[ I \ f G V2L 0 \ f V \ 

\ lQa\= \ K 0 ATa/T . (33) 

\lQb j Vo 0 K J \ An/T J 


Consider conversion of heat Iga into work —IV at ATa 7 ^ 
0 and ATf, = 0. The figure of merit and power factor are 


ZT = 


2L^ 


2aS^T 


GK-2L‘^ K-2aS‘^T’ 


P = 2aS^ 


(34) 


respectively, with k = Klu/{AT). The above hgure of 
merit and power factor are greater than those of the lon¬ 
gitudinal and transverse thermoelectric effects which are 


Z,T = ZfT = 


aS^T 

:-aS‘^T' 


Pi = Pt= aS^. 


(35) 


IV. CALCULATION OF THERMOELECTRIC 
PERFORMANCE 




E,o/(kBT) 


Eio/(kBT) 


(c) 



Eio/(kBT) 


(d) Pi 



-6 -3 0 3 6 

Eio/(kBT) 


(e) 


Pt (f) 


P 


ni 



Eio/(kBT) 


Eio/(kBT) 


FIG. 4. (Color online) Figures of merit ZiT (a), ZtT (b), and 
ZmT (c), and power factors (in unit of 10“^ W K~^) 
Pi (d), Pt (e), and Pm (f) versus the QD energy Eeo and 
Ero for Aqd = 20 meV and Kp = 0.05 W K“^. The 
temperature is T = 400 K, and the sheet density of QDs in 
each polymer layer (thickness 20 nm) is 4 x 10^^ cm“^. QD 
tunneling linewidth is F =30 meV. 

We now calculate the transport coefficients using ma¬ 
terial parameters. The robustness of the device per¬ 
formance is tested by including the randomness of QD 
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energy. Beside the transport mechanism described in 
Sec. II, there are other mechanisms conducting heat 
among reservoirs. The most important one is the heat 
conduction across the resisting layers by phonons. Poly¬ 
mers are good thermal insulators with heat conductivity 
0.05 ~ I W m“^ K“^.— Phonon thermal conductiv¬ 
ity should be much reduced in the nano-scale thin films 
concerned here due to abundant scattering with the em¬ 
bedded QDs and interfaces^ We take the phonon ther¬ 
mal conductivity as Kp = 0.05 W m“^ K~^ (e.g., bulk 
rayon^^ has thermal conductivity of 0.05 W m“^ K“^). 
The thermal conductance across a resisting layer is Kp = 
AupT/lqd. Adding this contribution to the transport 
equation and using Eq. ca) leads to 

Ki = ^+Kti + Kp, Ki 2 = - Kti - Kp, 

r2 

K2 = ^+Ktr + Kti + 2Kp. (36) 

Cr 

The variance of the QDs energy is several tens of meV as 
revealed by experiments.— The thickness of the resisting 
layer is taken as Iqd = 20 nm. The thickness of the source, 
the drain, and the cavity are all equal to 80 nm. High 
QD density, 2 x 10^^ cm“^, has been realized in polymer 
matrices of thickness around 20 nmj^ It is favorable to 
have more than one layer of QDs in each resisting layer 
to enable sufficient electron tunneling. Serial tunneling 
through several QDs may happen in the transmission 
across the resisting layers with many QDs.— We suggest 
to incorporate ~ 4 x 10^^ cm“^ QDs into a 20 nm poly¬ 
mer layer which corresponds to a volume density of QDs 
as ~ 2 X 10^® cm“® (average inter-dot distance ~ 8 nm). 
Taking into account of finite QD size, the inter-dot tun¬ 
neling linewidth is I ^ 30 meV depending on materials 
and structures which is taken as 30 meV here (see ex¬ 
perimental measured value of 30 meV in Ref. |4^ . And 
we take P = = rg(£;) = T^^E) = Td^jiE) = 

30 meV {\/i,j,E) as the QD tunneling linewidth used in 
calculating the tunneling rate. We shall study how the 
performance of the device varies with the variance and 
the mean value of the energy of QDs. The random QD 
energy in the left (right) resisting layer is modeled by a 
Gaussian distribution centered at E^o (A^o) with a vari¬ 
ance 
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with P = l,r. The QD energy and size can be con¬ 
trolled by various chemical^^ and physical^! methods 
during growth. We consider situations with temperature 
T = 400 K. The energy zero is set to be the equilibrium 
chemical potential. The band edge of the semiconductor 
that constitutes the source, drain, and cavity layers is 
200 meV below the chemical potential (a typical value 
for heavily-doped semiconductors). The electrical con¬ 
ductivity, thermopowers, and thermal conductivities are 


calculated according to Eqs. (El), (El), and El- Based 
on those transport coefficients we calculate the figures 
of merit and power factors, ZiT, ZtT, ZmT [Figs. HJa), 
SKb), andSKc)], Pi, Pt, and P^ [Figs. Hid),HKe), andSKf)]. 

Fig. m indicates that the figure of merit ZiT is opti¬ 
mized at the two points Eeo — ±3fc_BT with Ero = 0, 
while ZtT is optimized at Ego = —Ero — E^ksT. These 
results, which are consistent with the results in Ref. 113, 
can be understood as the balance between large electrical 
conductivity, large thermopower, and small thermal con¬ 
ductivity in optimizing the figure of merit. The power 
factors. Pi and Pt, are optimized at parameters simi¬ 
lar to those of ZiT and ZtT, respectively. We find that 
when ZtT is optimized, ZmT is only slightly larger than 
ZfT. For other situations, ZmT is considerably greater 
than both ZiT and ZtT. Particularly near the two points 
Eio = Ero ~ ESksT, the enhancement of figure of merit 
induced by cooperative effect is significant. For the power 
factor, cooperative effect always leads to considerable en¬ 
hancement of power factor, unless when Pi or Pt are close 
to zero. 

The above results reveal that cooperative effects can 
effectively improve the figure of merit and power factor 
for thermoelectrics in three-terminal geometry. Such im¬ 
provement is especially useful for systems of which the 
electronic structure has not been fully optimized. Hence 
cooperative effects offer an additional way to improve the 
performance of thermoelectrics that are potentially useful 
for realistic systems. We also note that the largest figure 
of merit and power factor for the longitudinal thermoelec¬ 
tric effect are ZiT = 0.6 and Pi = 1.6 x 10“® W m“^ K~^ 
respectively, while the largest figure of merit and power 
factor for the transverse thermoelectric effect is ZtT = 1 
and Pt = 4: X 10“® W m“^ K~^ respectively. This result 
confirms the conclusion in Refs. [l3, 123 - 1251 and that 
the transverse thermoelectric effect in three-terminal ge¬ 
ometry is of potential advantages. 

In order to check the robustness of the effect in real¬ 
istic situations we discuss the effect of the broadening of 
QD energy Aqd and the energy difference Ero — Ego for 
Ego = 110 meV (~ iksT). Increase of the broadening of 
QD energy reduces the thermopower and increases 
the electronic heat conductivity. Therefore the figures of 
merit and power factors for the longitudinal, transverse, 
and cooperative thermoelectric effects are all reduced. 
From Fig. O one finds that considerably large figures of 
merit and power factors can still be obtained for broad¬ 
ening of QD energy up to 50 meV. In experiments the 
full width at half-maximum of photoluminescence spectra 
can be as small as 35 meV (i.e., the variance is 15 meV) — 
Thus the proposed device is of potential application val¬ 
ues. Finally in the above discussions only one energy 
level in each QD is considered. Careful calculation with 
higher energy levels (100 meV higher) included indicates 
that the figure of merit and the power factor are even 
larger [see Appendix]. 

We also plot the figures of merit and power factors for 
longitudinal, transverse, and cooperative thermoelectric 










FIG. 5. (Color online) Figures of merit ZiT (a), ZtT (b), and 
ZmT (c), and power factors (in unit of 10~^ W K~^) Pi 
(d), Pt (e), and P^ (f) versus the variance of the QDs energy 
Aqd and the difference of the mean values of QD energy in 
the two resisting layers Ero — Eeo. Here Ero is varying while 
Eeo ~ 110 meV (~ SksT) is fixed. The temperature is T — 
400 K, Kp = 0.05 W m~^ K“^, and the sheet density of QDs 
in each polymer layer (of thickness 20 nm) is 4 x 10^^ cm“^. 
QD tunneling linewidth is F = 30 meV. 


FIG. 6. (Color online) Figures of merit ZiT (a), ZtT (b), 
and ZmT (c), and power factors (in unit of 10“"^ W m“^ K“^) 
Pi (d), Pt (e), and Pm (f) versus the tunneling linewidth of 
QDs F and the difference of the mean values of QD energy in 
the two resisting layers Ero — Eeo. Here Ero is varying while 
Eeo ~ 110 meV (~ SksT) is fixed. The temperature is T = 
400 K, Up = 0.05 W m~^ K“^, and the sheet density of QDs 
in each polymer layer (of thickness 20 nm) is 4 x 10^^ cm“^. 
The variance of the QDs energy is Agd = 20 meV. 


effects as functions of the tunneling linewidth of QDs F 
and the energy difference E^o — Eeo for Eeo = 110 meV 
(~ SksT) in Fig. [51 Unlike the monotonic dependence on 
the variance of the QDs energy Aqd, the hgures of merit 
and power factors first increases and then decreases with 
increasing F. This behavior is because at small F increase 
of F enhances electron tunneling and hence improves the 
electrical conductivity and the power factors. The en¬ 
hancement of electron tunneling also improves electronic 
heat conductivity and reduces the effect of phonon heat 
conductivity on the figures of merit. Therefore, the fig¬ 
ures of merit of the longitudinal, transverse, and cooper¬ 
ative thermoelectric effects are improved ad well. How¬ 
ever, the tunneling linewidth F also induces broadening 
of the energy of transported electron. When such broad¬ 
ening is comparable with or larger than the thermal en¬ 
ergy ksT, it considerably reduces the thermopowers and 


increases the electronic thermal conductivity, hence the 
power factors Pi, Pt, and Pm as well as the figures of 
merit ZiT, ZtT, and ZmT are reduced. The optimal tun¬ 
neling linewidths for ZiT, ZtT, and ZmT are F = 17.3, 
21.4, 21.4 meV respectively where the optimal figures of 
merit are ZiT = 0.668, ZtT = 0.987, and ZmT = 0.989 
respectively. Besides, the optimal tunneling linewidths 
for Pi, Pt, and Pm are F = 34.5, 39.4, and 38.6 meV 
respectively where the optimal power factors are Pi = 
1.43 X 10-3 W m-i K-2, Pt = 3.96 x IQ-^ W m'^ K-^, 
and Pm = 5.02 x IQ-^ W m-^ respectively. The op¬ 
timal hgure of merit and power factor for the cooperative 
thermoelectric effect are larger than those if the trans¬ 
verse and longitudinal thermoelectric effects. This result 
is consistent with the proof in Sec. m that the hgure 
of merit and the power factor of the cooperative ther¬ 
moelectric effect is greater than or equal to those of the 
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longitudinal and transverse thermoelectric effects. Mean¬ 
while the optimal performance of the transverse thermo¬ 
electric effect is also better than that of the longitudinal 
thermoelectric effect. Overall there are more parameter 
regions for the cooperative thermoelectric effect to have 
large values of the figure of merit and the power factor. 


(a) Zf^T / max(Z]T, ZfT) 
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FIG. 7. (Color online) The enhancement of figure of 
merit ZMT/max{ZiT, ZtT) versus the QD energy Eeo and 
Ero for Kp = 0.05 W m“^ and Agd = 20 meV (a), 
as well as when the phonon heat conductivity is increased 
to Kp = 1 W m~^ or when the QD energy broaden¬ 
ing is increased to Aqd = 150 meV (c). In (d) figures of 
merit ZmT as a function of Ero with Eeo = Ero is plotted 
for Kp = 0.05 W m“^ and Agd ~ 20 meV (curve with 
triangles), and for when the phonon heat conductivity is in¬ 
creased to Kp = 1 W m~^ (curve with dots) or when the 
QD energy broadening is increased to Agd = 150 meV (curve 
with squares). The temperature is T = 400 K, and the sheet 
density of QDs in each polymer layer (thickness 20 nm) is 
4 X 10^^ cm“^. QD tunneling linewidth is F = 30 meV. 


Finally we demonstrate the robustness of the coop¬ 
erative effects by examining the enhancement factor 
ZmT jTussiiZiT^ ZtT) of thermoelectric figure of merit for 
Kp = 0.05 W m“^ K~^ and = 20 meV, as well as 
when the phonon parasitic heat conductivity is increased 
to 1 W m“^ or when the QD energy broadening is 
increased to = 150 meV. The results are plotted for 
different QD energies in Fig. [3 Considerable enhance¬ 
ment of figure of merit by cooperative effect is found 
around Eeo ^ for a large portion of parameter re¬ 
gion. Moreover, this enhancement by cooperative effect is 
still effective when the phonon heat conductivity is much 
enhanced or when the QD energy broadening is signifi¬ 
cantly increased. This result reveals that the cooperative 
effect remains effective in improving thermoelectric effi¬ 
ciency even in systems with small figure of merit induced 
by significant parasitic heat conductivity [see Fig. [71[d)]. 


V. CONCLUSION AND DISCUSSIONS 

In summary we propose to enhance the thermoelectric 
efficiency and power by exploiting cooperative effects in 
three-terminal geometry. The three terminal geometry 
enables one electric current to couple with two temper¬ 
ature gradients with the help of inelastic transport pro¬ 
cesses. A scheme exploiting quantum-dots embedded in 
polymer matrices in multiple-layered structures is sug¬ 
gested to realize the principle. According to calculations 
based on material parameters, the figure of merit and 
power factor of the proposed structure are high, which 
indicates that layered resonant tunneling structures are 
potentially good thermoelectric systems i^^i^^ Marked im¬ 
provements of figure of merit and power factor by the 
cooperative thermoelectric effect are obtained. Remark¬ 
ably the enhancement of figure of merit and power factor 
induced by cooperative effects is robust to the parasitic 
phonon heat conductivity as well as quantum dots en¬ 
ergy broadening. Hence we shown that cooperative ef¬ 
fect offers an effective way to improve the figure of merit 
and power factor for three-terminal thermoelectric sys¬ 
tems, particularly useful for systems of which the elec¬ 
tronic structure has not been optimized. Study in this 
work indicates that exploiting geometric aspect, inelastic 
processes, and cooperative thermoelectric effects could 
provide alternative routes to high performance thermo¬ 
electrics. 
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APPENDIX: EFFECTS OF HIGHER LEVELS IN 
QDS ON THERMOELECTRIC PROPERTIES 

We calculate the figures of merit ZiT, ZtT, and ZmT 
and power factors Pi, Pt, and Pm for the situation when 
there is another energy level of 100 meV higher than the 
original one in each QD for both the left and the right 
polymer layers. The results are plotted in Fig. |8l We 
assume that the tunneling rate of the higher level is the 
same as the lower one, i.e., F = 30 meV. It is seen that 
the figures of merit as well as the power factors are all 
larger when the higher level is taken into account. Qual¬ 
itatively, the results here is similar to those in Fig. |T] but 
with the center shifted toward lower energy for both Eeo 
and Ero- This observation reveals that the main effect 
of the higher level is to enhance the electrical conduc¬ 
tivity as well as the thermopower, which can be under¬ 
stood easily since a higher energy channel is introduced. 
However, introducing such a channel also increases the 
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FIG. 8. (Color online) Figures of merit ZiT (a), ZtT (b), and 
ZmT (c), and power factors (in unit of W K“^) 

Pi (d), Pt (e), and Pm (f) versus the QD energy Eeo and Ero 
when there is another energy level of 100 meV higher than 
the original one and of the same tunneling rate in each QD 
for both the left and the right polymer layers. Aqd = 20 meV 
and Kp — 0.05 W m“^ K“^. The temperature is T = 400 K, 
and the sheet density of QDs in each polymer layer (thickness 
20 nm) is 4 X 10^^ cm“^. QD tunneling linewidth is F = 
30 meV. 


thermal conductivity which normally would reduce the 
figure of merit. 

To clarify the underlying mechanism, we plot the ther¬ 
mal conductivity ki as a function of Epq when E£q = ErQ. 
In Fig. |9] we plot three different contributions of ki: 


and 


bp . 


The 


aSfT, [a/ie^T)] [ ((i? - - {E - 

figure of merit 2(1 is given by 2(1 = aSlT /{ki — oSiT). 
Indeed the thermal conductivity due to energy uncer¬ 
tainty [cr/(e^r)] {^{E — ^ —{E — increases when 

higher level is introduced. However, in Fig.[9Ka) (with the 
higher level) when aSfT reaches its maximum value, it is 
very close to the other two contributions. In comparison, 
in Fig.lHKb) (without the higher level) when aSfT reaches 
its maximum value, it is considerably smaller than the 
phonon heat conductivity Hp. Therefore, the figure of 
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FIG. 9. (Color online) Heat conductivity rei (black full curve) 
and its three contributions: aS^T (green dashed curve), 

[cr/(e^r)] 1^ ((S — — {E — /i}^ j (blue dotted curve), and 

Kp (red chained cnrve) as functions of Ero = Eeo for (a) with 
the higher level and (b) without the higher level. The higher 
level is 100 meV higher than the original level in each QD. 
Aqd ~ 20 meV and Kp = 0.05 W m“^ K“^. The temperature 
is T = 400 K, and the sheet density of QDs in each poly¬ 
mer layer (thickness 20 nm) is 4 x 10^^ cm“^. QD tunneling 
linewidth is F = 30 meV. 


merit 2iT is enhanced when the higher level is taken 
into account. This particular feature is because in the 
case without the higher level, phonon heat conductivity 
predominately limits the figure of merit, rather than the 
variance of electronic energy. 


* jianhua.jiang.phys@gmail.com 

^ T. C. Harman and J. M. Honig, Thermoelectric and ther- 
momagnetic effects and applications (McGraw-Hill, New- 
York, 1967); H. J. Goldsmid, Introduction to Thermoelec¬ 
tricity (Springer, Heidelberg, 2009). 

^ G. J. Snyder and E. S. Toberer, Nat. Mater. 7 , 105 (2008); 
A. Shakouri, Ann. Rev. of Mater. Res. 41 , 399 (2011); T. 
M. Tritt, Ann. Rev. of Mater. Res. 41 , 433 (2011). 

® G. J. Snyder and T. S. Ursell, Phys. Rev. Lett. 91, 148301 
(2003). 

^ R. Venkatasubramanian, Phys. Rev. B 61, 3091 (2000); J.- 
K. Yu, S. Mitrovic, D. Tham, J. Varghese, and J. R. Heath, 


Nat. Nanotechnol. 5, 718 (2010); N. Nakpathomkun, H. 

Q. Xu, and H. Linke, Phys. Rev. B 82, 235428 (2010); 

R. Venkatasubramanian, E. Siivola, T. Golpitts, and B. 
O’Quinn, Nature 413 , 597 (2001). 

® A. I. Boukai et ah, Nature 451 , 168 (2008); B. Poudel et 
ah, Science 320 , 634 (2008); P. Pichanusakorn and P. Ban- 
daru. Mater. Sci. Eng. R-Rep. 67, 19 (2010); A. J. Min- 
nich, M. S. Dresselhaus, Z. F. Ren, and G. Chen, Energy 
Environ. Sci. 2, 466 (2009); C. J. Vineis, A. Shakouri, A. 
Majumdar, and M. G. Kanatzidis, Adv. Mater. 22, 3970 
(2010); Z.-G. Chen, G. Han, L. Yang, L. Gheng, and J. 
Zou, Prog. Nat, Prog. Nat. Sci. 22 , 535 (2012); J.-F. Li, 




































11 


W.-S. Liu, L.-D. Zhao, and M. Zhou, NPG Asia Mater. 2, 
152 (2010). 

® M. S. Dresselhaus, G. Chen, M. Y. Tang, R. Yang, H. 
Lee, D. Wang, Z. Ren, J.-P. Fleurial, and P. Gogna, Adv. 
Mater. 19, 1043 (2007). 

^ U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 

® L. D. Hicks and M. S. Dresselhaus, Phys. Rev. B 47, 12727 
(1993); ibid., 47, 16631 (1993). 

® G. D. Mahan and J. O. Solo, Proc. Natl. Acad. Sci. (USA) 
93, 7436 (1996). 

T. E. Humphrey, R. Newbury, R. P. Taylor, and H. Linke, 
Phys. Rev. Lett. 89, 116801 (2002); T. E. Humphrey and 
H. Linke, Phys. Rev. Lett. 94, 096601 (2005). 

P. Kim, L. Shi, A. Majumdar, and P. L. McEuen, Phys. 
Rev. Lett. 87, 215502 (2001); D. Li, Y. Wu, P. Kim, L. Shi, 
P. Yang, and A. Majumdar, Appl. Phys. Lett. 83, 2934 
(2003); J. Hone, M. Whitney, C. Piskoti, and A. Zettl, 
Phys. Rev. B 59, R2514 (1999); C. Yu, L. Shi, Z. Yao, D. 
Li, and A. Majumdar, Nano Lett. 5, 1842 (2005). 

F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, 
and J. P. Pekola, Rev. Mod. Phys. 78, 217 (2006). 
Thermoelectric Nanomaterials, edited by K. Koumoto and 
T. Mori (Springer-Verlag, Berlin, 2013). 

S. L. Korylyuk et al, Sov. Phys. Semicond. 7, 502 (1973); 
V. P. Babin et al, Sov. Phys. Semicond. 8, 478 (1974); A. 
Kyarad and H. Lengfellner, Appl. Phys. Lett. 89, 192103 
(2006); C. Reitmaier, F. Walther, and H. Lengfellner, 
Appl. Phys. A 99, 717 (2010); H. J. Goldsmid, J. Elec¬ 
tron. Mater. 40, 1254 (2011). 

S. Zippilli, G. Morigi, and A. Bachtold, Phys. Rev. Lett. 
102, 096804 (2009); B. Rutten, M. Esposito, and B. 
Cleuren, Phys. Rev. B 80, 235122 (2009). 

O. Entin-Wohlman, Y. Imry, and A. Aharony, Phys. Rev. 
B 82, 115314 (2010). 

R. Sanchez and M. Biittiker, Phys. Rev. B 83, 085428 
(2011); B. Sothmann, R. Sanchez, and A. N. Jordan, 
arXiv: 1406.5329 

B. Sothmann, R. Sanchez, A. N. Jordan, and M. Biittiker, 
Phys. Rev. B 85, 205301 (2012). 

J.-H. Jiang, O. Entin-Wohlman, and Y. Imry, Phys. Rev. 
B 85, 075412 (2012). 

T. Ruokola and T. Ojanen, Phys. Rev. B 86, 035454 

( 2012 ). 

B. Sothmann and M. Biittiker, Europhys. Lett. 99, 27001 

( 2012 ). 

L. Simine and D. Segal, Phys. Chem. Chem. Phys. 14, 
13820 (2012). 

J.-H. Jiang, O. Entin-Wohlman, and Y. Imry, Phys. Rev. 
B 87, 205420 (2013). 

A. N. Jordan, B. Sothmann, R. Sanchez, and M. Biittiker, 
Phys. Rev. B 87, 075312 (2013) 

J.-H. Jiang, O. Entin-Wohlman, and Y. Imry, New J. Phys. 
15, 075021 (2013). 

Y. Imry, O. Entin-Wohlman, and J. H. Jiang, WO Patent 
2,013,035,100. 


S. Juergens, F. Haupt, M. Moskalets, and J. Splettstoesser, 
Phys. Rev. B 87, 245423 (2013); C. Bergenfeldt, P. 
Samuelsson, B. Sothmann, C. Flindt, and M. Biittiker, 
Phys. Rev. Lett. 112, 076803 (2014); R. Bosisio, C. Gorini, 

G. Fleury, and J.-L. Pichard, New J. Phys. 16, 095005 
(2014); ibid., arXiv:1407.7020. 

B. Sothmann, R. Sanchez, A. N. Jordan, and M. Biittiker, 
New J. Phys. 15, 095021 (2013). 

H. L. Edwards, Q. Niu, and A. L. de Lozanne, Appl. Phys. 
Lett. 63, 1815 (1993); H. L. Edwards, Q. Niu, G. A. Geor- 
gakis, and A. L. de Lozanne, Phys. Rev. B 52, 5714 (1995). 
J. R. Prance, C. G. Smith, J. P. Griffiths, S. J. Chorley, 
D. Anderson, G. A. C. Jones, 1. Farrer, and D. A. Ritchie, 
Phys. Rev. Lett. 102, 146602 (2009). 

S.-K. Park, J. Tatebayashi, and Y. Arakawa, Appl. Phys. 
Lett. 84, 1877 (2004); S. Tonomura and K. Yamaguchi, 

J. Appl. Phys. 104, 054909 (2008); M. Jo, T. Mano, Y. 
Sakuma, and K. Sakoda, Appl. Phys. Lett. 100, 212113 
( 2012 ). 

Y. Yang, chap. 10 in Physical Properties of Polymers 
Handbook, edited by J. E. Mark (AlP Press, Woodbury, 
NY, 1996). 

W. R. Algar and U. J. Krull, Langmuir 24, 5514 (2008); 
S. Ishii, R. Ueji, S. Nakanishi, Y. Yoshida, H. Nagata, T. 
Itoh, M. Ishikawa, and V. Biju, J. Photochem. Photobiol. 
A 183, 285 (2006). 

K. Zhang, Y. Zhang, and S. Wang, Scientific Reports 3, 
3448 (2013). 

The polymer electronic band gap should be considerably 
larger than the semiconductor band gap, so that the charge 
carriers in the doped semiconductor regions (i.e., the cav¬ 
ity, source, and drain) will not diffuse into the polymer 
resisting layers. 

R. S. Whitney, Phys. Rev. Lett. 112, 130601 (2014). 

M. Gong, K. Duan, C.-F. Li, R. Magri, G. A. Narvaez, and 

L. He, Phys. Rev. B 77, 045326 (2008); A. D. Lad and S. 
Mahamuni, Phys. Rev. B 78, 125421 (2008). 

H. Haug and A. P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors, (Springer, Berlin, 1996). 

L. Onsager, Phys. Rev. 37, 405 (1931); 38, 2265 (1931). 
Y. Demirel and S. 1. Sandler, J. Phys. Chem. B 108, 31 
(2004). 

H. T. Odum and R. C. Pinkerton, Am. Sci. 43, 331 (1955). 
J.-H. Jiang, Phys. Rev. E 90, 042126 (2014). 

C. Van den Broeck, Phys. Rev. Lett. 95, 190602 (2005); 

M. Esposito, K. Lindenberg, and C. Van den Broeck, Phys. 
Rev. Lett. 102, 130602 (2009). 

S. A. Gurvitz, Phys. Rev. B 57, 6602 (1998). 

J. Berezovsky, O. Gywat, F. Meier, D. Battaglia, X. Peng, 
and D. D. Awschalom, Nat. Phys. 2, 831 (2006). 

K. Koumoto et al, in CRC thermoelectric handbook: Micro 
to Nano, D. M, Rowe ed. (CRC press, USA, 2006). 

A. Agarwal and B. Muralidharan, Appl. Phys. Lett. 105, 
013104 (2014). 


